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1  Introduction 


In  the  1940’s,  Kovasznay  [1]  was  interested  in  obtaining  measurements  of  fluid  turbulence 
quantities  in  the  wake  of  artillery  projectiles.  His  technique  used  duplicate  slides  of  a  spark 
range  shadowgraph  (see  Section  3)  and  a  light  source  to  perform  a  spatial  correlation  of 
the  light  intensity  in  the  wake  imaged  on  the  shadowgraph.  He  then  used  the  theory  of 
Taylor  [2]  to  obtain  statistical  properties  of  the  three  dimensional  density  fluctuations  in  the 
wake  from  the  two  dimensional  light  intensity  fluctuations  in  the  shadowgraph.  This  was  a 
slow  process.  With  the  use  of  digital  scanning  devices  and  computers,  the  process  is  much 
quicker.  This  report  shows  the  results  of  using  modern  equipment  interrogate  a  spark 
range  shadowgraph  with  Kovasznay’s  technique. 

2  Background 

Ballistic  ranges  and  wind  tunnels  have  been  used  for  several  decades  to  determine  the  force 
and  moment  coefficients  that  describe  the  stability  and  the  motion  of  a  projectile  in  free  flight. 
For  many  years,  these  facilities  provided  the  only  means  of  determining  those  coefficients.  In 
the  past  ten  years,  computational  fluid  dynamics  has  developed  enough  to  allow  the  value  of 
these  coefficients  to  be  determined  from  a  computer  solution  to  the  Navier-Stokes  equations. 
Nevertheless,  the  computer  solutions  need  information  about  the  viscous  layer  of  the  flow 
field  surrounding  the  projectile,  especially  in  transonic  flight. 

Computer  solutions  to  the  Navier-Stokes  equations  generally  use  one  of  two  closure 
schemes:  an  algebraic  model  to  define  a  turbulent  viscosity  or  a  model  of  the  turbulent 
kinetic  energy  production  and  dissipation.  Both  models  benefit  from  experimental  flow  field 
data.  For  example,  the  “Baldwin-Lomax”  model  [3]  is  used  successfully  in  computational 
fluid  dynamics  for  closure  to  the  Navier-Stokes  equations.  The  Baldwin-Lomax  model  takes 
advantage  of  the  Prandtl-Van  Driest  formulation  for  the  turbulent  viscosity.  In  this  model, 
the  mixing  length,  varies  with  distance  from  a  boundary  (or  wall)  according  to: 

^  = /cy[l -exp(-y+/yl+)].  (1) 

In  this  equation,  «  and  are  constants,  y  is  the  dimensional  distance  from  the  boundary, 
and  y"*"  is  a  non-dimensional  distance  from  the  boundary.  The  mixing  length  translates  to 
a  physical  length  within  the  flow  field,  namely  the  size  of  the  eddies  that  are  transporting 
the  bulk  of  the  momentum.  Experimental  measurements  of  the  size  of  these  eddies  would 
benefit  the  modeling.  Similarly,  measurements  of  the  turbulent  energy  spectra  in  the  flow 
field  help  the  production-dissipation  models  of  turbulence. 

The  computer  solutions  calculate  the  velocity  and  pressure  fields  in  the  boundary  layer 
surrounding  a  projectile  in  flight.  Once  the  pressure  at  the  surface  of  the  projectile  is 
calculated,  the  forces  and  moments  (and,  finally,  the  coefficients  used  in  trajectory  and 
stability  analysis)  can  be  calci  lated.  While  in  a  wind  tunnel,  techniques  such  as  hot  wire 
anemometry  or  laser-Doppler  velocimetry  measure  velocities  in  the  flow  field,  the  same  is 
generally  not  true  for  ballistic  ranges.  Similarly,  pressure  measurements  in  the  flow  field  can 
be  made  in  a  wind  tunnel,  but  are  not  made  routinely  in  ballistic  ranges.  Instead,  ballistic 
ranges  are  used  to  find  the  position  of  a  projectile  at  a  certain  time  and  then,  with  the 
equations  of  motion,  that  information  is  used  to  calculate  the  ballistic  coefficients. 
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Detailed  experimental  information  about  the  flow  field  around  a  body  is  usually  only 
available  from  wind  tunnel  measurements,  and  not  from  ballistic  ranges.  However,  wind 
tunnel  measurements  of  velocity  and  pressure  fields  are  not  ideal,  since  the  projectile  is 
not  in  free  flight  and  the  Reynolds  number  can  be  significantly  different  from  free  flight 
conditions.  Wind  tunnel  measurements  can  also  be  subject  to  unknown  disturbances  from 
the  walls  and  mounting  structures.  Thus,  any  technique  for  extracting  flow  field  information 
from  the  ballistic  range  would  be  of  use  to  those  who  model  the  flow  field  on  the  computer. 

The  point-to-point  statistical  record  of  the  velocity  fluctuations  in  a  turbulent  flow  field 
available  from  laser-Doppler  velocimetry  is  nearly  impossible  to  obtain  in  a  ballistic  range; 
even  if  one  measurement  could  be  made,  the  number  of  firings  required  for  an  ensemble 
average  would  be  prohibitive.  Large  field  optical  techniques  would  be  of  great  interest  if  they 
were  able  to  capture  the  instantaneous  two  dimensional  flow  field  around  a  projectile.  In 
this  way,  an  ensemble  average  might  be  made  in  ten  or  so  firings.  The  technique  discussed 
in  this  paper  is  very  far  from  being  able  to  give  a  quantitative  picture  of  the  flow  field 
surrounding  the  projectile.  Nevertheless,  the  use  of  simple  shadowgraphs  and  computers 
suggests  that  more  information  can  be  squeezed  from  the  ballistic  range  that  is  of  importance 
to  computational  modeling. 

3  The  Shadowgraph 

Shadowgraphs  can  be  produced  in  the  Transonic  Range  at  the  Ballistic  Research  Laboratory. 
It  is  a  207  meter  long  building  with  a  6  meter  square  test  section  instrumented  with  25 
orthogonal  shadowgraph  stations,  as  described  by  Rogers  [4].  In  addition  to  the  orthogonal 
shadowgraph  stations,  it  is  possible  to  obtain  a  direct  photograph  of  the  flow  field  surrounding 
the  projectile  through  the  use  of  a  mosaic  station.  This  station  consists  of  a  platform 
positioned  near  the  line  of  flight  of  the  projectile  and  covered  with  sheets  of  photographic  film. 
(The  mosaic  is  a  direct  shadowgraph;  the  orthogonal  stations  produce  indirect  shadowgraphs, 
i.e.  the  shadowgraph  image  is  projected  onto  a  screen  and  the  screen  is  photographed.  See 
Rogers  for  the  details.)  When  the  projectile  is  above  the  platform,  a  light  source  is  triggered 
to  obtain  the  shadowgraph.  This  composite  picture  generally  magnifies  the  flow  field.  The 
mosaic  station  for  the  present  work  wcis  placed  at  a  distance  of  122  meters  from  the  gun; 
Figure  1  is  a  schematic  of  the  station. 

Figure  2  shows  a  shadowgraph  obtained  from  the  mosaic  station;  a  schematic  is  presented 
in  Figure  3  to  clarify  what  is  shown  in  the  shadowgraph.  The  projectile  in  the  shadowgraph 
is  an  155mm,  M864  Base-Burn  projectile;  as  shown  in  the  schematic,  only  a  small  portion 
of  the  projectile  base  is  imaged.  There  is  some  gas  injection  shown  in  the  shadowgraph, 
although  the  propellant  in  the  M864  has  not  reached  a  steady-state  burn. 

Shadow  patterns  on  the  shadowgraph  are  produced  by  a  redistribution  of  the  intensity 
of  the  illuminating  spark  light  source.  The  rays  of  light  are  deflected  as  though  there  were 
small  lenses  in  the  wake  region;  hence,  the  formation  of  a  “shadow”  depends  on  the  changes 
in  deflection,  not  on  the  actual  value  of  the  deflection  angle.  For  small  angles  of  deflection. 

Ax  =  W  (2) 

where  the  displacement  of  the  light  ray  on  the  shadowgraph  is  Ax,  /  is  the  distance  between 
the  centerline  of  the  wake  and  the  shadowgraph,  and  0  is  the  angle  of  deflection  at  the  wake. 
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Figure  1;  Shadowgraph  Arrangement  in  the  Transonic  Range 


Figure  2;  Schematic  of  Shadowgraph  of  M864  Base-Burn  Projectile  with  Gas  Injection 
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Uniform  deflection  would  move  the  rays  of  light  uniformly;  only  changes  in  6  will  result 
in  a  varying  Ax  at  the  shadowgraph.  These  changes  in  Ax  correspond  to  local  vaxiations 
of  the  light  intensity  on  the  screen,  and  thus  appear  as  lighter  and  darker  regions  on  the 
shadowgraph. 

Density  variations  in  the  wake  account  for  the  effective  lenses  since,  from  the  empirical 
law  due  to  Gladstone  and  Dale  [5], 


-  =  K  or  n  =  pK  +  l  (3) 

P 


where  K  is  the  Gladstone-Dale  constant,  6  is  equal  to  n  —  1  (n  is  the  refractive  index*)  and 
p  is  the  density  of  the  gas.  For  air,  the  density  varies  with  the  temperature,  T,  and  the 
pressure,  p,  according  to: 

Tq  p  f  . 

P  =  (4) 

1  Po 


where  the  zero  subscript  refers  to  an  initial  value.  Hence,  variations  of  the  pressure  or 
temperature  in  the  wake  of  the  projectile  result  in  density  variations  in  the  wake,  and, 
thus,  refractive  index  variations.  This  maJces  shadows  on  a  surface  opposite  the  light  source 
illuminating  the  wake. 

If  a  density  fluctuation  creates  an  apparent  convex  lens,  then  there  is  the  possibility  that 
the  light  rays  passing  through  the  apparent  lens  will  focus  before  they  reach  the  film.  In 
this  case  the  resulting  shadowgraph  will  not  have  a  one-to-one  correspondence  between  its 
shadows  and  the  density  fluctuations.  The  shortest  focal  length  that  could  result  from  the 
shadowgraph  used  in  this  study  can  be  determined  from  the  lens  maker’s  equation: 


]_ 

f 


(5) 


where  r,-  are  the  radii  of  curvature  of  each  side  of  the  lens  and  /  is  the  focal  length  (positive 
if  on  the  real  side  of  the  lens,  negative  if  on  the  virtual  or  image  side).  The  worst  case  would 
occur  with  a  flat  surface  on  the  incoming  ray  side  of  the  lens  and  a  radius  of  the  resolution 
of  the  film  (300  lines  per  millimeter)  on  the  outgoing  side.  Equations  3  and  4  can  be  used 
to  determine  the  change  in  the  refractive  index  of  air  given  that  the  pressure  is  twenty  per 
cent  lower  in  the  wake  than  in  the  free  stream  (a  worst  Ccise).  With  these  conditions,  the 
focal  length  is  about  2  meters.  The  wake  was  about  one  meter  away  from  the  film  in  this 
study;  therefore,  it  is  assumed  that  a  one-to-one  correspondence  exists  between  the  density 
fluctuations  and  the  record  on  the  shadowgraph. 

A  simple  relation  holds  for  the  deflection  angle  of  the  light  rays  with  changes  in  refractive 
index,  namely: 

J-z  dx 
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-dz 


(6) 


where  z  is  the  direction  normal  to  the  shadowgraph  and  the  wake  is  2Z  thick.  (The  coordinate 
system  is  shown  in  Figure  4.)  Thus,  summing  all  the  variations  in  refractive  index  through 
the  thickness  of  the  wake  yields  the  resulting  angle  of  deflection  of  the  incoming  ray  of  light. 

‘The  refractive  index  is  the  ratio  of  the  velocity  of  light  in  vacuum  to  the  velocity  of  light  in  the  given  medium.  The  ratio 
of  the  difference  in  velocities  to  the  velocity  in  the  medium  is  the  value  of  £. 
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Figure  3;  Shadowgraph  and  Schematic  of  M864  Base-Burn  Projectile 
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Figure  4:  Distance  and  Sign  Conventions 
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The  first  step  in  analyzing  the  shadowgraph  is  to  determine  the  relationship  between  the 
intensity  of  the  light  recorded  on  film  (the  illumination)  and  the  density  fluctuations  in  the 
waJce.  The  cheinge  in  illumination,  /,  (energy /area/time)  is  in  proportion  to  the  rate  of 
chajige  of  the  deflection  of  light  rays  with  distance  [6]: 


AI  de 

T 


(7) 


This  is  true  as  long  as  the  refractive  index  gradients  are  not  too  large,  as  is  the  case  for 
shadowgraphs  of  the  waJje  of  projectiles.  The  same  is  not  necessarily  true  for  flames.  This  can 
have  consequences  for  powered  projectiles,  such  as  rocket-assisted  or  base-burn  projectiles. 
Combining  Equations  3,  6,  and  7  and  taking  three  dimensions  into  consideration  yields: 


AI 

I 


IK  f+z  d^p  d^p  ,  IK 
Po  J-z  ax^  ay^  po 


(8) 


where  I  remains  the  distance  from  the  centerline  of  the  wake  to  the  shadowgraph,  K  is  the 
Gladstone- Dale  constant,  po  is  the  air  density  in  the  free  stream,  p  is  p{x,y,z)  and  y  is  the 
direction  in  the  same  plane  as  and  perpendicular  to  x.  At  this  point  it  is  convenient  to  define 
a  dimensionless  parameter  that  indicates  the  illumination  on  the  shadowgraph: 


F{x,y) 


Iq  -  /(x,  y) 

i{x,y) 


AI 

I 


(9) 


where  the  zero  subscript  refers  to  the  undisturbed  illumination  of  the  shadowgraph.  Function 
F  can  be  negative;  it  is  not  a  gray  scale. 


4  Spatial  Correlation 


The  correlation  coefficient  of  the  density  fluctuations  in  the  wake  of  a  projectile  was  deter¬ 
mined  by  Kovasznay  [1]  using  an  optical  technique.  He  made  a  duplicate  of  the  shadowgraph 
transparency,  placed  the  duplicate  behind  the  original,  and  used  a  traversing  mechanism  to 
move  one  shadowgraph  with  respect  to  the  other  in  the  plane  of  the  shadowgraphs.  Us¬ 
ing  a  light  source  and  a  photodetector,  he  could  measure  the  transmittance  through  both 
shadowgraphs  resulting  from  different  separations  (the  distance  between  a  corner  of  one  shad¬ 
owgraph  and  the  same  corner  of  the  other).  This  formed  the  local  correlation  coefficient,  J, 
from  a  separation,  s: 

Jt{s)  =  J  hr,, T{x)T{x  +  s)dx  (10) 

where  hnt  is  the  intensity  of  the  light  source  and  x  is  the  vector  locating  position  in  the 
overlapping  areas  of  the  shadowgraphs.  Term  T  is  the  transmissivity  of  the  shadowgraph, 
which  is  related  to  F{x,y).  This  correlation  coefficient  for  the  two  dimensional  shadowgraph 
was  then  related  to  the  three  dimensional  density  fluctuations  by  Equation  8: 


Jt{s)(xJ  p\^pdz 


(11) 
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The  local  correlation  of  the  density  fluctuations  is  defined  as: 

J,{s)  =  J  p{X)p{X  +  s)dX  (12) 

where  X  is  a  three-dimensional  position  vector,  so  that  Equation  1 1  can  be  rewritten  as: 

Jr(s)  =  CJ^  [V^J  J,{s)  dz'  dz  (13) 

where  C  is  a  constajit.  Since  Jp  is  a  function  of  the  displacement  s  only,  and  p  is  statistically 
homogeneous  [7],  the  Laplacian  is  now: 


= 


iA 

sds 


(14) 


By  defining  a  new  variable  ^  =  z  —  z',  Equation  13  becomes: 


'  (2Z-0  Jp{s)d^ 
0  ^  ^ 


(15) 


If  the  thickness  of  the  wake  is  much  larger  than  the  largest  turbulent  scales  that  can  be 
resolved  in  the  correlation,  then  the  length  2Z  is  effectively  infinite.  In  this  case.  Equation  15 
caji  be  solved  by  inverse  Fourier  transform. 


5  Discrete  Interrogation  of  the  Shadowgraph 

Another  technique  for  spatial  correlation  of  the  shadowgraph  is  to  discretize  the  illumination 
with  the  use  of  an  optical  scanning  device.  In  this  operation,  the  shadowgraph  is  divided 
into  small  square  areeis,  (Ajj),  «is  shown  in  Figure  5.  The  illumination  in  a  portion  of  each 
area  is  averaged;  the  value  of  the  average  is  represented  by  7,^  .  Figure  5  depicts  the  portion 
used  in  the  averaging  process.^  An  area’s  average  in  this  work  could  take  on  integer  values 
between  zero  (for  no  illumination)  and  255  (for  100%  illumination). 

Local  correlation  of  the  illumination  can  be  represented  by: 


Ma,b) 


J+b 


(16) 


where  a  and  b  are  the  number  of  row  and  column  separations.  Hence,  the  duplicate  shadow¬ 
graph  technique  described  above  can  be  replicated  by  a  “shift  and  multiply”  technique  on 
the  computer  which,  if  done  in  Fourier  space,  is  very  fast. 

The  density  field  in  the  wake  can  be  transformed  using  the  following  expansion: 


P{x,y,z)=  JL 


|p|<f 


kl<f  H<f 


^  px  qy  rz' 

a:  +  i:+2z, 


(17) 


In  this  equation,  P,  Q,  and  R  are  the  number  of  Fourier  modes  in  the  x,  y,  and  z  directions, 
function  g  is  the  spectrum  of  the  density,  and  G(a)  =  cos(q)  -f  sin(a).  Terms  A*  and  Ay  are 

^ li,]  !■  obtained  from  a  subset  of  Ajj  due  to  the  hardware  being  used  (see  next  section  for  details). 
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Figure  5:  Schematic  of  Discrete  Values  of  the  Illumination  on  Shadowgraphs 


the  lengths  of  the  wake  as  imaged  on  the  shadowgraph  in  the  subscripted  direction.  When 
this  expansion  is  substituted  into  Equation  8,  the  factor  ~~  results  from  the  differentia¬ 
tion.  In  the  limit  of  Z  ^  oo  (which  is  reasonable  for  the  shadowgraph  where  the  scale  of 
the  turbulence  patterns  is  much  smaller  than  the  thickness  of  the  wake)  and  treating  the 
summation  over  r  as  infinite,  then: 


and  the  final  result  is: 

F{x,y) 


=  TT 


TT  /  ic 

|p|<|  ki<# 


(18) 


(19) 


6  Equipment 

The  digital  scanning  device  used  was  an  Eikonix  Model  1412  Digital  Imaging  Camera,  built 
by  the  Eikonix  Imaging  Systems  division  of  Kodak.  Since  this  was  the  first  project  in  the 
BRL  to  attempt  to  make  quantitative  measurements  with  this  device,  considerable  effort 
was  put  into  the  calibration  of  the  equipment  and  the  determination  of  error  thresholds  for 
the  system.  Several  errors  and  important  omissions  were  found  in  the  documentation  from 
Eikonix;  hence,  the  following  is  presented  as  an  aid  to  future  users  of  this  device. 

The  camera  consists  of  a  55mm  f/2.8  Nikon  lens  in  front  of  a  4096  element  charged 
coupled  device  (CCD)  array.  This  one-dimensional  CCD  array  can  be  stepped  through  4096 
positions  allowing  up  to  a  4096x4096  pixel  scan.  The  elements  of  this  array  are  five  microns 
wide  by  seven  microns  high,  with  a  two  micron  gap  between  each  element  in  the  five  micron 
direction.  Each  sample,  or  pixel,  is  quantized  by  a  12-bit  analog  to  digital  (A/D)  converter. 
The  finite  size  of  each  pixel  means  that  the  output  of  the  scanner  is  a  point  sampling  of  the 
input  image  after  a  rectangular  “box  car”  filter  has  been  applied. 
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Figure  6:  Detail  of  Uncalibrated  Scanner  Output 


The  object  being  scanned  can  either  be  transmissive  (e.g.  photographic  negatives)  or 
opaque  (e.g.  prints).  A  fluorescent  light  table  is  used  for  the  former  case,  while  four  tungsten 
lamps  provide  top  lighting  for  opaque  objects.  The  exposure  is  controlled  by  both  the  f-stop 
of  the  lens  and  the  “integration  time”  of  the  CCD  array,  which  is  discussed  below.  A 
programmable  Alter  wheel  allows  for  color  scanning  via  red,  green,  and  blue  filters;  but,  only 
gray  scale  was  used  and  discussed  in  this  report. 

Since  the  spacing  of  the  CCD  array  is  fixed,  the  resolution  of  a  scan  is  controlled  by  the 
distance  of  the  camera  above  the  scanned  surface.  At  its  closest  focus,  the  camera  yields 
approximately  a  55mm  wide  scanning  area  over  the  4096  element  array,  or  a  13.4  micron 
sample  spacing  (1892  saimples  per  inch).  Higher  resolution  could  be  obtained  by  changing 
lenses  on  the  camera  (Eikonix  offers  an  optional  105mm  lens).  A  practical  limit  for  the 
largest  single  area  that  can  be  scanned  is  approximately  1  meter  side  to  side,  which  yields  a 
sample  spacing  of  244  microns  (104  samples  per  inch).  For  the  55mm  lens,  the  width  of  the 
scanned  area  was  found  to  be  related  to  camera  position  by: 

scan  width  =  4.95  camera  position  -f  13mm  (20) 

The  camera  is  connected  to,  and  controlled  by,  a  Sun  Microsystems  3/280  computer  via  an 
IEEE-488  bus  interface.  The  light  sources,  camera  distance,  focus,  and  f-stop  are  controlled 
manually.  Scanning  area,  exposure,  CCD  array  normalization,  etc.  are  computer  controlled. 
Data  can  be  transferred  to  the  computer  as  either  8  or  12  bits  per  sample  (in  the  former 
caise,  the  lower  4  bits  of  the  A/D  converter  are  ignored). 


Naive  use  of  the  scanner  can  result  in  fairly  substantial  errors.  Figure  6  shows  a  contrast 
enhanced  result  of  a  scan  of  a  uniform  white  field.  The  brightness  variations  have  been 
exaggerated  by  a  factor  of  fifteen.  The  lines  running  vertically  through  the  image  result 
from  the  non-uniform  response  of  the  CCD  array  elements.  The  horizontal  lines  in  the 
image  result  from  brightness  variations  in  the  light  source  during  the  scanning  process.  The 
following  discussion  shows  that  corrective  measures  can  be  taken  to  reduce,  but  not  eliminate, 
both  of  these  errors. 

Each  element  of  the  CCD  array  acts  as  an  integrator  of  the  incoming  light.  Between  the 
time  that  an  element  is  cleared  or  “set  to  zero”  and  the  time  that  it  is  sampled  by  the  A/D 
converter  it  will  sum  up  the  illumination  that  falls  on  it.  This  time  interval  is  called  the 
“integration  time.”  Each  element  of  the  array  responds  linearly  to  the  incident  illumination. 
However,  the  constants  of  this  linear  equation  are  different  for  every  element,  and  will  vary 
with  temperature. 

Ideally,  if  an  element  received  no  light  during  its  integration  time  the  sampled  output 
would  read  zero.  In  practice,  each  element  will  yield  some  small  non-zero  value  referred  to 
as  the  “dark  current.”  If  exposed  to  equal  amounts  of  light,  the  sampled  outputs  of  every 
element  should  read  the  same;  but,  due  to  sensitivity  variations  from  element  to  element, 
this  does  not  happen.  To  “normalize”  the  response  of  the  array,  a  bias  and  gain  value  should 
be  computed  for  every  element  of  the  array  and  then  scanned  values  should  be  corrected  by 
these  factors. 

The  Eikonix  camera  contains  bias  and  gain  buffers  which  can  be  loaded  with  these  cor¬ 
rection  factors.  This  caJibration  procedure  involves  exposing  the  array  to  black  (e.g.  with 
the  lens  cap  on)  to  find  bias  values,  and  then  exposing  the  array  to  a  uniform  white  field 
to  find  gain  values.  To  achieve  the  best  possible  uniform  white,  the  authors  suggest  fully 
defocusing  the  camera  over  either  the  cleaned  surface  of  the  light  table  (for  transmissive 
scans)  or  a  plain  white  sheet  of  paper  (for  reflective  scans).  Observation  of  a  pin  hole  light 
source  showed  the  point  spread  when  fully  defocused  to  be  approximately  one  twentieth  of 
the  entire  scanning  area. 

Variation  in  the  light  source  during  integration  can  also  be  a  problem  and  was  the  source 
of  the  horizontal  lines  in  Figure  6.  The  major  source  of  this  variation  is  the  120Hz  “flicker” 
of  the  lights  caused  by  the  AC  power  supply.  It  is  important  that  the  CCD  array  be  exposed 
for  an  integreil  multiple  of  this  flicker  time,  otherwise  the  amount  of  light  received  during 
the  integration  time  will  depend  on  the  phase  of  the  lights  at  the  instant  each  scan  line 
exposure  begins.  The  manual  makes  the  erroneous  claim  that  the  flicker  of  the  tungsten 
lights  is  somehow  taken  into  account  by  the  camera.  While  the  rather  massive  filaments  and 
the  low  operating  voltage  of  the  tungsten  lamps  reduces  this  problem,  flicker  is  still  present. 
A  photocell  was  used  to  measure  this  light  variation  and  was  found  to  be  slightly  over  two 
per  cent  peak  to  peak. 

The  flicker  of  the  fluorescent  light  table  is  extreme,  since  the  light  almost  completely  shuts 
off  between  cycles.  Getting  the  correct  integration  time  is  critical.^  Given  the  potential 
aisymmetries  in  the  fluorescent  discharge  cycle,  the  flicker  cycle  rate  should  be  treated  as 


^There  are  several  incorrect  or  confusing  formulas  for  the  integration  time  in  the  manual.  For  the  model  1412  camera  it  is 
simply 

clock  rate  count  •  integration  time  count 
tntegratton  time  =  - g -  microseconds 
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60Hz  (16.67  msec).  There  was  one  other  low  frequency  error  associated  with  the  fluorescent 
light  table  that  could  be  seen  with  the  naJced  eye  and  is  impossible  to  remove;  this  error 
results  from  unstable  oscillation  of  the  gas  plasma. 

One  additional  error  associated  with  the  Eikonix  camera  was  found.  After  calibrating  the 
bias  and  gain  tables,  individual  scan  lines  still  show  a  systematic  offset  between  gray  values 
reported  for  even  and  odd  pixels.  The  CCD  array  used  in  the  camera  is  designed  where 
even  and  odd  elements  axe  separated  into  two  distinct  output  stages.  It  was  not  determined 
whether  this  error  resulted  from  a  flaw  in  the  automatic  caJibration  sequence,  or  from  an 
instability  in  the  CCD  device.  The  end  result  was  a  small  noise  component  at  the  Nyquist 
rate  (wave  number  of  one  half  sample)  in  the  horizontal  direction. 

For  the  scans  in  this  report,  prints  of  the  shadowgraph  were  used  with  the  reflective 
tungsten  lighting.  An  integration  time  of  33.333  msec  was  used.  From  one  scan  to  the  next, 
sample  values  varied  by  about  five  per  cent.  This  seems  disappointingly  high,  but  is  of  the 
same  order  as  the  quantization  noise  for  our  images. 

7  Results 

The  shadowgraph  of  Figure  3  was  digitized;  the  result  is  shown  in  Figure  7.  This  figure  is 
from  the  lower  left  of  the  four  plates  in  the  shadowgraph  configuration.  The  dark  area  at  the 
upper  right  of  the  shadowgraph  is  the  base  of  the  projectile.  Two  areas  were  selected  from 
this  shadowgraph  for  statistical  evaluation  and  are  shown  by  the  boxes  in  the  figure.  Area 
one  is  the  “course  grained”  region  near  the  base  of  the  projectile,  and  area  two  is  in  the  finer 
grained  region  to  the  lower  left.  Each  of  these  areas  were  digitized  with  nearly  the  maximum 
resolution  of  the  scanner  and  contain  4096  by  4096  samples.  Each  sample  is  14.6  microns 
on  the  photographic  print.  Careful  measurement  of  the  projectile  in  the  figure  showed  that 
the  shadowgraph  was  magnified  by  a  factor  of  1.63.  The  projectile  wcis  traveling  at  a  Mach 
number  of  about  1.2. 

Detailed  pictures  of  these  two  areais  are  shown  in  Figures  8  and  10.  Each  of  these  pictures 
covers  «in  area  of  36.2  mm  on  a  side;  each  sample  was  thus  8.84  microns  on  a  side  (113 
pixels/mm).  Statistical  quantities,  showing  the  distribution  of  gray  values  among  pixels,  are 
shown  in  Table  1.  The  associated  histograms  are  shown  in  Figures  9  and  11.  The  statistical 
quantities  are  similar  in  both  areas,  although  area  one  is  skewed  to  a  greater  extent.  Both 
areas  use  less  than  half  of  the  gray  scale  resolution  available.  This  would  suggest  that  an 
additional!  6  dB  of  signal  to  quantization  noise  ratio  (SNR)  might  be  obtained  by  either  a 
higher  contrast  print,  or  a  better  scanner  exposure. 

Figures  12  through  15  show  the  results  from  shift  and  multiply  in  one  row  or  column  of 
the  digitized  areas.  These  are  normalized  such  that  the  maximum  is  unity  in  the  unshifted 
position.  The  mean  value  for  the  entire  area  has  been  subtracted.  At  larger  values  of  shift, 
less  samples  are  included  in  the  autocorrelation  estimate.  As  a  result,  the  variance  of  the 
estimate  goes  up.  To  avoid  misinterpretation  of  this  data,  the  biased  autocorrelation  estimate 
is  shown. 

After  the  initial  decrease  from  the  unshifted  position,  all  four  correlations  display  a  sinu¬ 
soidal  variation  for  about  20  mm.  Periodicity  in  these  results  is  an  indication  of  the  size  of 
the  scales  of  motion  in  the  wake.  These  figures  show  that  the  most  signiflcajit  scale  that  is 
resolved  from  this  process  has  a  length  of  about  one  millimeter;  this  is  the  same  length  that 
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Kovasznay[l]  found  in  his  work.  This  is  based  on  the  first  zero  crossing  in  each  of  the  figures, 
since  this  is  the  distance  at  which  patterns  on  the  shadowgraph  are  uncorrelated.  Because 
the  autocorrelation  was  performed  using  a  zero-based  average  for  the  gray  value,  negative 
values  of  the  autocorrelation  are  possible  and  indicate  the  distance  at  which  patterns  with 
gray  vadues  below  the  average  are  correlated  with  those  above.  Area  two  has  a  somewhat 
higher  frequency  of  the  variations  than  area  one,  and  this  is  consistent  with  the  notion  of  a 
finer  greiin  in  the  turbulent  structure  of  area  two. 

Figures  16  through  19  are  the  averages  of  the  autocorrelation  estimates  of  all  the  rows 
and  all  the  columns  in  the  two  areas.  The  autocorrelation  approaches  zero  near  a  separation 
distance  of  one  millimeter  in  all  four  figures.  In  the  average,  the  periodicity  that  was  so 
apparent  in  the  single  row  or  column  autocorrelation  has  been  smoothed.  Now  the  autocor¬ 
relation  appears  to  be  that  for  a  random  function:  rapid  decay  from  a  maximum  value  to 
zero.  This  shows  that  while  in  any  given  row  or  column  there  are  particular  frequencies  that 
are  dominant,  those  frequencies  change  a  short  distance  away.  Averaged  over  the  entire  area, 
a  broad  rainge  of  frequencies  (or  length  scales)  are  present  in  a  continuous  manner.  This  is 
consistent  with  the  dynamics  of  turbulent  wake  flow. 

The  results  of  the  two-dimensional  autocorrelations  are  shown  in  Figures  20  through  23; 
both  a  contour  plot  and  a  gray  scale  image  are  presented.  The  values  of  the  autocorrelation 
are  indicated  on  the  contour  plot;  the  gray  scale  image  shows  positive  values  as  levels  of  gray, 
with  black  indicating  values  of  zero  or  less.  These  were  computed  over  an  area  of  512  square 
pixels;  the  shift  was  512  pixels  in  both  the  column  and  row  directions.  The  512  pixels  chosen 
were  from  the  lower  left  corner  of  each  area.  The  local  minimum  in  area  one  occurs  about 
1.5  mm  from  zero  separation;  it  is  skewed  such  that  it  is  farther  from  zero  in  the  column 
direction.  In  area  two,  the  skew  is  toward  the  row  direction;  however,  the  magnitude  of  the 
distance  is  once  again  near  1.5  mm.  Area  two  appears  to  have  some  coherence  along  the 
diagonal  (see  Figure  10)  and  this  shows  up  in  the  two-dimensional  correlation  as  a  diagonal 
streak  of  white  <ifter  the  local  minimum. 

Rotational  correlations  give  an  indication  of  length  scale  in  the  wake,  and  these  are  shown 
with  five  degrees  of  rotation  in  Figures  24  and  25.  Rings  form  about  three  millimeters  from 
the  center  of  the  picture,  and  these  rings  have  a  spacing  of  about  one  millimeter.  Again,  this 
suggests  the  correlation  of  density  fluctuations  in  a  very  graphic  manner. 

The  energy  spectra  of  the  two  areas  (using  both  row-  and  column-averaged  shift  and  mul¬ 
tiply)  tire  shown  in  Figure  26.  This  spectra  is  a  correlogram,  which  is  the  Fourier  transform 
of  the  autocorrelation  estimates.  Several  references  [2,  8,  9,  10]  show  estimations  of  what  the 
spectra  should  look  like.  Typically  on  such  plots  the  energy  is  normalized  by  the  product 
of  the  fluid  viscosity  (to  the  |)  and  the  dissipation  (to  the  quarter  power),  and  the  wave 
number  is  multiplied  by  an  appropriate  length  scale  (either  the  Taylor  micro-scale  or  the 
Kolmogorov  length  scale).  The  result  is  a  curve  with  a  slope  of  -1  out  to  moderate  wave 
numbers,  and  then  a  change  to  a  slope  of  —  |  for  large  wave  numbers  for  wall  bounded  shear 
flows.  In  wakes,  the  slopes  are  expected  to  be  similar.  In  an  isotropic  turbulence,  the  —  | 
slope  dominates  the  wavenumber  space  [2].  This  is  shown  schematically  in  Figure  27.  The 
change  in  slope  occurs  at  unity  in  both  the  dimensionless  wavenumber  and  dimensionless 
energy.  (The  initial  portion  of  the  curve  is  nearly  flat.) 

In  the  present  work,  the  values  are  not  normalized,  since  it  was  unclear  which  length  scale 
to  use.  None-the-less,  on  a  log-log  scale  both  normalizations  produce  an  offset  only,  so  the 
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slope  of  the  curve  will  remain  the  same  no  matter  which  scale  is  used.  The  results  show 
two  distinct  linear  regions,  with  the  slope  for  higher  wave  numbers  more  negative  than  at 
low  wave  numbers.  At  low  wave  numbers  the  slope  for  both  digitized  areas  is  roughly  -1,  in 
good  agreement  with  predicted  values.  The  above  is  true  for  wave  numbers  up  to  0.2;  at  this 
point,  there  is  another  change  in  slope.  The  data  become  somewhat  more  noisy,  also.  These 
data  appear  to  be  accurate,  since  the  large  spike  at  the  highest  wave  number  corresponds 
to  the  Nyquist  rate  error  in  the  scanner  CCD  array  described  in  the  Equipment  section.  At 
present,  there  is  no  explanation  for  the  second  change  in  slope.  There  is  some  offset  between 
the  two  digitized  areas,  and  this  suggests  (as  shown  above)  different  length  scales  in  the  two 
areas. 
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Figure  7:  Digitized  Shadowgraph,  Areas  Studied 


Table  1:  Statistics  for  the  Digitized  Shadowgraphs 


Quantity! 

Area  One 

Area  Two 

Minimum  g 

14 

16 

Maximum  g 

140 

137 

Number  of  g 

127 

122 

Equivalent  dB 

44 

43 

Mean  g 

48 

58 

Variance 

143 

155 

Standard  Deviation 

15.6 

12.4 

Skewness 

4615 

1579 

1  fTerm  g  is  the  gray  scale  value 

14 


Per  cent 


Figure  8:  Area  One,  Expanded  View 
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Figure  9;  Histogram  of  Area  One 
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Per  cent 


Figure  10:  Area  Two,  Expanded  View 
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'igure  16:  1-D  Autocorrelation,  Row  Average  in  Area  One 


;ure  17:  1-D  Autocorrelation,  Column  Average  in  Area  One 
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Y  Shift  (mm) 


Figure  20:  2-D  Autocorrelation  Contours,  Area  One 


Figure  21:  2-D  Autocorrelation  in  Gray  Scale,  Area  One 
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Figure  22:  2-D  Autocorrelation  Contours,  Area  Two 


Figure  23:  2-D  Autocorrelation  in  Gray  Scale,  Area  Two 
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Figure  24;  Circular  Correlation  of  Five  Degrees,  Area  One 


Figure  25:  Circular  Correlation  of  Five  Degrees,  Area  Two 
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Figure  26;  Energy  Spectra  of  the  Row  and  Column  Averaged  Data  in  Areas  1  and  2 
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Figure  27:  Schematic  of  the  Energy  Spectra  for  Wall  Bounded  Shear  Flow 
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Conclusions  and  Future  Work 


The  conclusions  reached  in  this  work  axe  as  follows: 

•  Discrete  interrogation  of  shadowgraphs  gives  results  that  are  nearly  identical  to  analog 
techniques. 

•  Discrete  interrogation  is  much  faster  than  analog  techniques  and  it  allows  the  data  to 
be  manipulated  easily. 

•  Further  investigation  should  begin  with  a  homogeneous,  turbulent  wake  with  known 
maximum  length  scales. 

•  Mathematical  relations  between  length  scales  and  the  density  fluctuations  must  be 
established. 

The  goal  of  this  work  was  to  help  computational  fluid  mechanics  by  gaining  knowledge 
about  the  scales  of  motion  in  a  turbulent  flow  field.  This  first  effort  is  far  short  of  that  goal, 
but  establishes  a  collection  of  tools  that  may  help  to  exploit  ballistic  range  shadowgraphs  for 
this  purpose.  There  axe  three  efforts  that  would  enhance  the  technique  without  developing 
complicated  three-dimensional  image  equipment. 

First,  a  conditioned  averaging  approach  to  the  simple  technique  described  in  this  report 
would  give  a  better  indication  of  the  length  scale.  The  conditions  would  pick  the  size  of  the 
area  to  be  averaged  such  that  the  random  nature  of  the  flow  field  would  not  dominate  the 
results  and  the  small  scale  coherence  could  be  traced. 

Second,  in  an  effort  to  study  the  far  field  wake  or  the  wake  at  different  angles  of  yaw, 
several  mosaic  stations  could  be  set  up  to  image  the  wake.  The  far  field  wake  could  be  imaged 
by  triggering  three  or  four  mosaic  stations  at  the  same  time,  thus  capturing  a  significant 
portion  of  the  wake  in  three  or  four  locations.  Triggering  the  stations  as  the  projectile 
passed  through  would  show  the  near  field  wake  under  different  conditions  and,  using  several 
cissumptions  and  Taylor’s  theory,  might  allow  approximations  of  three-dimensional  structure. 

Finally,  a  two  slice  mosaic  station  could  be  constructed  to  edlow  two  shadowgraphs  to  be 
obtained  at  right  angles.  This  might  yield  some  information  about  the  three-dimensional 
structure  of  the  wake. 
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Computer 

Compile  Options 

Shifts/sec 

MFLOPS 

X  VAX 

VAX  11/780 

cc  -0 

.42 

0.13 

1.0 

Sun  3/280 

cc  -0  -f68881 

.42 

0.13 

1.0 

Sun  3/280 

gcc  -0 

.70 

0.22 

1.7 

Gould  9080 

cc  -0 

1.35 

0.43 

3.2 

VAX  6320 

cc  -0 

1.71 

0.55 

4.1 

Alliamt  FX/8 

cc  -0 

1.72 

0.55 

4.1 

SGI  4D/240 

cc  -03 

8.90 

2.80 

21.0 

Cray  X-MP/48 

cc 

13.20 

4.20 

31.0 

IBM  R6000 

cc  -0 

27.00 

8.60 

64.0 

Cray  2 

see 

189.00 

60.00 

450.0 

Cray  X-MP/48 

cc  -holeveU  -hivdep 

233.00 

75.00 

554.0 

Cray  X-MP/48 

see 

345.00 

110.00 

821.0 

Table  2:  Two  Dimensional  Autocorrelation  Computation  Speed 


Appendix  A.  Computational  Speed 

Most  of  the  computations  performed  in  this  report  were  of  a  very  simple  nature;  however,  the 
volume  of  data  involved  demamds  substantial  processing  power.  While  a  detailed  comparison 
of  different  computers  used  in  the  laboratory  is  beyond  the  scope  of  this  report,  we  present 
one  example  of  the  relative  speeds  of  various  machines  for  the  typical  kind  of  processing  this 
report  required. 

The  two  dimensional  autocorrelation  involves  point  multiplying  two  copies  of  an  image, 
one  of  which  is  shifted  in  the  horizontal  and/or  vertical  direction.  The  results  of  these 
multiplies  are  then  summed  together.  A  simple  C  program  to  compute  this  autocorrelation 
was  written,  the  inner  loop  of  which  is  essentially  a  vector  dot  product.  In  the  report,  512 
shifts  in  width  and  height  (262144  total)  were  used,  each  multiplied  and  summed  over  an 
cirea  of  512  by  512  samples.  The  total  number  of  floating  point  adds  and  multiplies  required 
is  137  billion.^ 

Table  2  compares  the  rates  at  which  a  selection  of  computers  could  perform  a  simple  two 
dimensional  autocorrelation.  The  type  of  computer,  along  with  the  compiler  and  options 
used  is  shown.  The  value  of  shifts  per  seconds  is  for  a  400  by  400  sample  area.  Millions 
of  floating  point  operations  per  second  (MFLOPS)  is  found  by  multiplying  the  shifts  per 
second  by  0.32,  as  320  thousand  multiplies  and  adds  are  required  for  each  shift.  Finally,  the 
speed  normalized  to  that  of  the  VAX  11/780  is  shown.  In  all  ca.ses  the  computation  was 
done  using  full  precision  (64  bit). 

A  few  characteristics  are  worth  pointing  out.  The  Sun  3/280  is  shown  as  speeding  up 
by  more  than  fifty  per  cent  by  the  change  in  C  compilers.  The  Cray  X-MP  also  shows  a 
near  fifty  per  cent  increase  under  different  compilers.  More  interesting  is  the  two  cases  of 
the  Cray  X-MP  using  cc.  With  the  option  to  “ignore  vector  dependencies”  the  inner  loop  of 
the  autocorrelation  code  vectorized  and  resulted  in  an  18  times  speedup.  This  was  typical 

*Ufe  of  the  Feet  Fourier  lyanifonn  to  perform  the  autocorreletion  could  reduce  the  operation  count  by  approximately  a 
factor  of  fifty;  but,  the  relative  requirements  discuued  here  still  apply. 
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of  several  programs  being  applied  to  this  data:  the  strictly  scalar  floating  point  speed  of  the 
X-MP.  was  about  30  times  a  VAX  11/780,  yet  vectorization  could  result  in  another  20  to  30 
fold  speedup. 

To  put  the  scale  of  the  above  numbers  in  better  perspective,  consider  the  original  two 
dimensional  autocorrelations  computed  in  this  report.  On  the  Sun  3/280  where  the  images 
were  scanned,  one  autocorrelation  would  have  taken  over  seven  days  of  computer  time.  On 
the  Cray  X-MP/48,  these  could  each  be  done  in  twenty  minutes. 
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Appendix  B.  Turbulent  Wake  Flow 

Wake  flows  show  similarity  in  patterns  at  different  sections  downstream  of  the  projectile. 
Self-preservation  of  these  patterns  moving  downstream  with  the  main  flow  is  because  the 
un-bounded  turbulent  region  in  the  wake  is  narrow  with  respect  to  its  length.  Streamwise 
speed  is  at  least  an  order  of  magnitude  greater  than  transverse  speed,  but  the  turbulence 
intensity  is  only  twice  as  great  in  the  streamwise  direction.  In  this  condition,  changes  in  flow 
patterns  axe  not  exaggerated  in  the  streamwise  direction  (hence,  they  are  “self-preserving”), 
but  axe  strongly  dependent  on  their  history  [2]. 

Self-preservation  would  only  hold  for  the  middle  range  of  wavenumbers  of  the  turbulence 
(see  Section  7)  without  production  of  turbulence;  a  grid-turbulence  in  a  wind  tunnel  is  an 
example  of  self-preservation  in  a  decaying  turbulent  flow  field.  With  turbulence  production, 
self-preservation  might  not  occur,  since  the  scales  of  motion  would  be  dependent  on  the 
turbulence  production.  In  a  wake,  however,  the  production  is  dependent  on  the  turbulence 
produced  upstream  and  transported  downstream  by  diffusion  and  convection  [8].  Thus,  self- 
preservation  might  be  expected  to  occur  and  there  is  strong  evidence  of  this  in  the  similarity 
of  turbulence  patterns  in  the  image  of  wake  flows. 

Similarity  in  a  flow  field’s  patterns  means  that  all  velocities  are  represented  by  one  char¬ 
acteristic  velocity  and  the  scales  of  motion  are  represented  by  a  characteristic  length.  Using 
the  similaxity  assumption,  which  is  strongly  dependent  on  self-preservation,  the  equations  of 
motion  for  a  turbulent  wake  can  be  simplified.  Unfortunately,  most  experimental  measure¬ 
ments  [2]  in  axi-symmetric  wake  flows  suggest  that  similarity  is  at  best  “incomplete,”  and 
this  is  true  only  in  the  fax  field  wake. 

Given: 

1.  the  difficulty  in  developing  adequate  turbulence  models  to  compute  the  wake  flow  field 
of  artillery  and  direct-fire  projectiles; 

2.  the  inability  to  simplify  the  equations  of  motion  in  the  near  wake; 

3.  the  lack  of  copious  experimental  data  (ais  in  wall-bounded  flows);  and 

4.  the  lack  of  development  of  experimental  techniques  that  could  be  used  to  investigate 
projectile  wake  flows 

it  is  very  important  to  develop  the  techniques  and  to  obtain  experimental  meaisurements  of 
turbulence  quantities  in  the  wake. 
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